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MULTIPOLE VORTEX BLOBS (MVB): 
SYMPLECTIC GEOMETRY AND DYNAMICS 

DARRYL D. HOLM & HENRY O. JACOBS 

DEPARTMENT OF MATHEMATICS 
IMPERIAL COLLEGE, LONDON SW7 2AZ, UK 


Abstract. Vortex blob methods are typically characterized by a regulariza¬ 
tion length scale, below which the dynamics are trivial for isolated blobs. In 
this article we observe that the dynamics need not be trivial if one is willing 
to consider distributional derivatives of Dirac delta functionals as valid vortic- 
ity distributions. More specifically, a new singular vortex theory is presented 
for regularized Euler fluid equations of ideal incompressible flow in the plane. 
We determine the conditions under which such regularized Euler fluid equations 
may admit vorticity singularities which are stronger than delta functions, e.g., 
derivatives of delta functions. We also describe the symplectic geometry associ¬ 
ated to these augmented vortex structures and we characterize the dynamics as 
Hamiltonian. Applications to the design of numerical methods similar to vortex 
blob methods are also discussed. Such findings illuminate the rich dynamics 
which occur below the regularization length scale and enlighten our perspective 
on the potential for regularized fluid models to capture multiscale phenomena. 
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1. Introduction 


Vortices are important in hydrodynamics because they are the sources for the 
incompressible flow held. The vorticity distribution at any instant of time deter¬ 
mines both the current state of the how and its future evolution, for given boundary 
conditions. This property holds for any Hamiltonian system, and it can indeed 
be shown that the dynamics of vortices can be usefully expressed in Hamiltonian 
form. In the vorticity and stream function formulation of an ideal incompressible 
planar huid, the evolution of the vorticity distribution u){x,y,t) is given by 

(1) dtU) — {u, -0} = dtu — dxoo dy^jj + dyU dx^J = 0 , 

where u = —A'ljj is the vorticity, ip is the stream function, and A = dxx + dyy is the 
Laplace operator. The corresponding (x, y) components of the Eulerian velocity 
held are given by 


{u,v) = {dyip, -dxip). 


If one is willing to view the vorticity a; as a distribution, one can consider point 
vortex solutions. In particular, point vortices are obtained if one considers the 
vorticity solution ansatz 



where Ti{t) E M., z = {x, y) G and 5z^{t) is the Dirac delta distribution centered 
at the point Zi{t) = {xi{t), yi{t)) e at a given time f G M. Substitution of this 
ansatz into Q yields the following well known hnite dimensional system in the 
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form of Hamilton’s canonical eqnations, 


( 2 ) 


^ = 0 , ip{z,t) = '^Ti{t)G{z - Zi{t)), 

i 

^ = dy'il){zi) , ^ = -d^tlj(zi) , 


where G{z) = ln(||z||) is the Green’s fnnction for the planar Laplacian. 

A point-vortex approximation to a continuons distribntion of vorticity for Euler’s 
fluid equations is problematic, though, because a point vortex induces a flow veloc¬ 
ity which becomes unbounded. However, when the point vortex is made smooth 
and bounded (regularized) the approximation becomes reasonable |Cho73] . 

For example, one may consider the regularized form of the vorticity equation 
given by choosing a translationally and rotationally invariant smoothing kernel 
of width 5 > 0 and dehning the regularized vorticity as Ks * u = —Aip while 
continuing to use ([^ to evolve ui in time. For example, Ks{z) = exp(—||2:|p/5^) is 
considered in |BM85j . In this case the point vortex ansatz yields (|^ again, except 
that the singular Green’s function G is replaced by the smooth kernel 


( 3 ) 


Gs{z) := Ks * G{z) = — (Ei(-||^||7^2) _ 2 ln(||^| 


where Ei(-) denotes the exponential integral function. The vorticity kernel Gs has 
no singularity at the origin for 5 > 0, and is known as a vortex blob. This system is 
the starting point for the vortex blob method, introduced in |Gho73] (albeit with 
a different regularization). 

The economy of the vortex blob method derives from the property that Dirac 
delta distributions are hyper-local (i.e. parametrized by position), and the prop¬ 
erty that the vorticity equation 0 admits Dirac delta distributions as solutions. 
However, there are many distributions which are localized to a similar degree (e.g. 
derivatives of delta functions, 9x7 )• 

In this paper, we study the more general vorticity solution ansatz. 


w(z, t) = 5 ^ rr‘(«)srs^L. 

i,m,n 


We hnd that this ansatz yields a closed hnite dimensional system which generalizes 
vortex blobs. We call these new carriers of vorticity multipole vortex blobs or 
MVBs. 


1.1. Main contributions. 

(1) Section briefly reviews the background for vortex methods in fluid mod¬ 
eling. 

(2) Section reviews the relationship between regularized fluids and vortex 
blob methods. 








4 


D.D. HOLM & H.O. JACOBS 


(3) Section derives the equations of motion for point vortices and MVBs as 
exact solutions of a regularized vorticity equation. 

(4) Section derives the conservation laws for these equations, such as energy, 
linear momentum, and angular momentum, and circulation. The derivation 
of these conserved quantities as symplectic momentum maps can be found 
in Appendix 

(5) Section explains the relationship between the dynamical systems for 
MVBs and an implicitly defined closed dynamical system which governs 
the spatial moments of the vorticity distribution. 

(6) Section discusses numerical aspects of using MVBs to model fluid dy¬ 
namics, such as approximations of initial conditions (subsection |7.2[ ), and 
grouping of computational nodes (subsection |7.1[ ). 

(7) Section presents the results of several numerical experiments involving 
small numbers of vortices, for 77 = 1,2, and 3. 

(8) MVB dynamics is Hamiltonian. We present the symplectic and Hamilton¬ 
ian structure of MVB dynamics in Section 

2. Background 

Vortex methods for fluid modeling predate the computer age and references to 
them can be found in the work of Helmholtz |LS111 see the introductory section] 
For example, the use of point vortices as idealized solutions can already be found 
in a 1931 paper concerning a “line of discontinuity” in planar fluid flow |Ros31] . At 
the beginning of their development, the inhnite velocities (and energies) associated 
to point vortices caused great difficulties, both numerically and theoretically. In 
fact, the point vortex approach did not produce a competitive numerical method 
until the 1970s, when the problems related to singularities were overcome by regu¬ 
larizing the singular vortex kernel to form a vortex blob. Stochastic perturbations 
were further included to model viscosity |Cho73] . These adjustments to the clas¬ 
sical point vortex method yielded the vortex blob method, which quickly became 
of practical use for realistic fluid flow modeling. In particular, the regularized 
system proved more amenable to error analysis. It was shown that the solutions 
of the vortex blob method converge to solutions of the Navier-Stokes equations 
in jHa179j . Later, stronger convergence rates were achieved by judicious choice of 
vortex kernels. By convolving the singular vortex kernel with sums of Gaussian 
smoothing kernels, a sequence of vortex blob kernels with faster convergence rates 
was found. Specifically, the convergence rate of the mth kernel was found to be of 
order for any q G (0,1) where h = 6'^ is a grid-spacing parameter and 5 > 0 is 
a length scale associated to the smoothing kernel |BM821 IBM85j . 

Simultaneously, the symplectic geometry of point vortices was clarihed in |MW83] 
by invoking Arnold’s interpretation of ideal fluids |Arn66] . The hndings of |MW83j 
were developed further in |GBV12j to handle fluid flow on manifolds with nontriv¬ 
ial homology. While this theoretical development clarihed the geometry of point 
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vortices, vortex blobs were sometimes thought to be a numerical “trick” which vi¬ 
olated the geometric interpretation. However, this thought was banished with the 
invention of the Euler-a model, a regularized model of ideal fluids with a parameter 
a representing the typical correlation length of fluctuations away from the mean of 
a Lagrangian fluid path |FHT01j . In particular, vortex blob solutions associated to 
a specihc kernel serve as exact solutions to the Euler-a model |OSni] . The Euler-a 
kernel is different from the kernels used in |Cho73] and |BM85] . A comparison of 
the Euler-a kernel to the m = 1 kernel of |BM85j is given in [RNPObj for vortex 
filament and vortex sheet motion. 

While vortex blobs performed well, they did not capture all of the qualitative 
richness observed in fluid vorticity dynamics. In particular, blobs of vorticity in 
real ideal fluids are known to change shape and deviate from initially circular dis¬ 
tributions. A numerical method is proposed in |Ros971 IRosOSj to capture these 
shape dynamics by adding basis functions with non-trivial moments in the study 
of vortex merger (see for example |MZM88l ILDVn2i IMEDLOb] ). Another distinct 
model obtained by projection onto a Hermite basis is described in [NSUW09]. This 
projection yielded a finite-dimensional systems which modeled the (truncated) mo¬ 
ments of the vorticity of an ideal incompressible fluid. The derivation of simplihed 
combinatorial formulas invoked by the dynamics of this model were discovered in 
[UWB12] and these formulas have made the method numerically tractable for a 
large number of moments. 

A dual approach to the moment based methods of the previous paragraph 
|Ros971 IRos05[ INSUW09] is to consider multipole based methods. This is the 
approach proposed in |Nic86] , where an initial vortex ansatz consisting of sums of 
distributional derivatives of dirac delta distributions is considered. Such an idea 
has occured intermittently in various forms in the literature, over many years. For 
example, a regularized vortex blob model, in the spirit of |BM82l IBM85] . which 
considered vorticity distributions of the form oj = + ^idx^zi + 

was investigated in |CN88j . Here it was proven that this augmentation of the 
traditional vortex method will yield faster spectral convergence than that of a tra¬ 
ditional vortex blob method. The current article considers higher order derivatives 
and can be seen as a natural follow up to |CN88] . More recently, dynamics have 
been derived for interactions of pure vortices and pure dipoles. These come from 
vorticity distributions of the form u = Yl,i + ^\dy5z^ ^) with 

the assumption that the locations of the dipoles and the vortices never overlap and 
that their self-interaction terms may be neglected |YTKn91 ITYKllj . In a differ¬ 
ent approach, approximations of dipoles are created by holonomically constraining 
vortices of opposite strength to be a hxed distance from one another, |TKN12] . 
The question remains, however, to what extent the dynamics of |TKN12] approx¬ 
imates those of |YTK09[ ITYKll] after self-interaction terms have been neglected. 
In summary, the removal of self-interaction terms is one of the primary obstacles 
to obtaining a multipole based generalization of the point vortex method |LSllj . 
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Moreover, the spectral error decay rates found in |Hal79[ IBM82t IBM85| ICN88j 
arise from the use of vortex blobs in place of (singular) point vortices. In this 
article we will follow this regularization based approach. 


3. Vortex blobs and regularized fluid models 

In this section we review a class of regularized fluid models and their relationship 
with vortex blob methods (for a more detailed discussion see |HNPn6j b The sort 
of fluid models we consider take the form 

dtu + u ■ Vu = 0 

oj = curl(LQ, ■ u). 

Where is a SE(2) invariant linear psuedo-differential operator with a length- 
scale parameter a > 0 such that limo-s.oQop = 1- When La is just the identity, 
the above “model” is Euler’s equations of motion for an ideal fluid. When La = 
(1 —where A is the Laplace operator, then we obtain the the Euler-a model, 
the solutions of which will converge to solutions of Euler’s ideal fluid equations as 
a > 0 vanishes IFHTOlj . 

We may replace u with its stream function, in order to rewrite the above 
equations as 

(4) dtu: + {V’, cu} = 0 

(5) (u = A(L„ • p)) 

This allows us to represent planar fluid dynamics in terms of scalar functions and 
distributions rather than vector-fields. 

The relationship between these regularized models and vortex blobs methods 
comes from first considering the point-vortex ansatz 

t) = ^Ti6{z - Zi{t)). 
i 

If the operator, Ao La, has a non-singular Green’s function, Ga, then substituting 
the ansatz into ([^ implies that 

(6) xpiz] t) = y^^TiGajz - Zi{t)) 

i 

We should note that when La is the identity (i.e. for an Euler fluid), then Ga is 
singular, and an extra argument (perhaps a physical one) must be presented in 
order to allow the resulting singular velocity helds. In this paper no such issues 
with singularity arise because we are modelling an Euler fluid with a regularized 
fluid where La has a non-singular Green’s function. 
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Substitution of ip into (|^ then implies the following equations of motion for the 
vortex cores Zi = {xi,yi) and the strengths Ti{t): 

(7) ^ = dyip(zi(t);t) , ^ =-d^ip(zi(t);t) , ^ = 0- 

When a = 0 and La = 1 this is nothing but the point-vortex method. When a > 0 
it is possible for ip to be much more regular, and we obtain various vortex-blob 
methods. In particular, we obtain the smooth vortex blobs of |BM85] . 

It is notable that @ and Q together form a finite dimensional ODE. The 
solutions of this ODE are exact solutions to the regularized fluid model. Again, 
this is valuable because the solutions of many regularized fluid models have been 
shown to converge to solutions of the ideal fluid equations as a vanishes. This 
paper seeks to generalize these point-like solutions to regularized fluid models to 
obtain a richer class of solutions with richer conservation properties. 


4. Equations oe motion 


In this section we derive the equations of motion for the time-dependent param¬ 
eters which specify multipole vortex blobs (MVBs). The zeroth order MVBs are 
just standard vortex blobs and the resulting equations of motion are those of the 
standard (non-stochastic) vortex blob algorithm |Cho73j . The hrst order MVBs 
are regularized dipoles and the equations of motion are those of |CN88] . Here we 
will derive the equations of motion for Nth order MVBs following the approach of 

[CN88]. 

Consider the ansatz for the vorticity, 

(8) = E rr(«)973p,, 

i^S m-\-n<N 


for spatially constant dynamical variables G M for i E S where S is some 

countable set. The stream function is 

(9) iiz,t) = Y. E - u«))■ 

i^S va-\-n<N 

The corresponding velocity field is given by 


( 10 ) 


U 


t) = t) = E, - z,{t )), 


v(z,t) = -d^iiz,t) = -Ei 


£S,m-\-n<N 


rr{t)d:p^+^d:;Gs{z - z,{t)). 


Examples of the types of velocity fields produced are depicted in figures [^through 
1^ on page 

We seek equations of motion for the r™”(t)’s and Zi{tys such that the velocity 
held (10) satishes the vorticity equation ([^. In the following calculations, we will 
not show the explicit time dependence of the dynamical variables. 
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Figure 1. A 0th order MVB with z = 0 and F = 1, using the 
kernel Gs of equation ([^. This form of the kernel produces one 
of the vortex blobs presented in |BM85] and the resulting numerical 
method is identical. 



Figure 2. The flow held around a 1st order MVB with T = 0, T^ = 
1, T^ = 1 is that of a regularized dipole. 



Figure 3. Various second order MVBs with all F’s set to zero ex¬ 
cept. (left) F^^ = 1, (middle) F^^ = 1, (right) F*^ = F^^ = 1. 
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We now find 
dtUJ= Y. 


ies 

m+n<N 


J'pmn 

dt 


p^m+l r^nx _ Tmn ^Vi ptmptn+1 x 

^zi j- i dt ® y ’ 


dt 


and 


dy'lljd^U= Y 


ies 

m+n<N 


By invoking 


of Appendix]^ we can rearrange the previous equation to obtain 
dy'iij d^uj = Y rr(-1)”^+”+'+'+" ^ 


ieS 
m+n<N 
i,k 


Similarly, we find 
d^tlj dyUj = ^ r™ 


■ 1 ) 


ies 

m+n<N 

l,k 


m+n+l+i+k 


n + 1 
k 




Substitution of these expressions into Q yields the vanishing of a linear com¬ 
bination of the distributions d^dl^Sz for m + n < A^ + 1. Since each of these 
distributions is linearly independent of the others (assuming the ZiS are distinct), 
their individual coefficients must each vanish independently. If we focus on the 
terms of the sum where m + n = A^ we obtain coefficients for and dyd at the 
core locations. The vanishing of these coefficients yields the dynamics for MVB 
cores 


( 11 ) 


= dyip{zi) , ^ = -d^i){zi). 


The vanishing of the coefficient of dz^ yields 

dr°’° 


dt 


= 0 . 


For ^ + k < N the vanishing of the coefficient of d^d'^Szi yields 


( 12 ) 


dFf 

dt 

(_l)£+fc 

E rr 


m>i 

n>k 

n-\-m<N 


m 

i-1 


+ 


n 

k-1 


m 


dm—l+l,n—k+l^ (^0 


We observe that dV^^/dt depends on 'ip at the vortex cores Zi, and the vortex core 
dynamics depend on ip as well. Fortunately, we already found that ip is purely a 
function of Fp and Zi, as stated in (|^. Thus ( p)^ and ( [ll| ) (with ([^) form a closed 
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finite dimensional system. Most notably, by construction the vorticity equation 
Q admits the A^th order MVB ansatz for the vorticity in ([^ as a solution when 
the Zi{t)^s and the rj(t)’s satisfy the just derived hnite dimensional system. 

Remark 4.1. In the point vortex method (i.e. the un-regularized case where = 
1), this derivation of the dynamics requires an extra step. In particular, one must 
discard the self-interaction term, which we will describe here. For point vortices, 
oj = —A'lp. Substituting the point-vortex ansatz u = ~ ^4^)) the 

equations of motion (|^ would then yield the non-sensical equation 

i* = V-^ Fj log \zi - Zj 

We say “non-sensical” because the right hand side explodes when you evaluate 
the ith term in the sum, the self-interaction term. Historically, it is customary 
to discard this self-interaction term based on physical and symmetry principles 
|MP941 Chapter 4]. In contrast, for blob methods the logarithmic kernel is replaced 
with a differentiable kernel function, such as a Gaussian. This allows one to 
retain the self-interaction terms. In the case of standard vortex blobs (i.e. 0th 
order MVBs), this distinction makes no difference because the gradient of the 
kernel vanishes at the origin and the self-interaction term contributes nothing to 
the dynamics. However, the derivatives of the kernel of degree 2 and higher do 
not vanish at the origin. As a result, the self-interaction terms do contribute to 
the dynamics for MVBs of order 2 and higher. The choice to discard the self¬ 
interaction terms in |YTK09j . versus our choice to keep them, explains one of 
the major discrepancies between our work and |YTK09] . In particular, |YTK09j 
was concerned with generalizing the (un-regularized) point vortex method in the 
same way that we have generalized the vortex-blob method. Once the ansatz 
u = ^ F™”9™ 9 ^ 5 ( 2 ; — Zi) was substituted into the equations of motion, they 
discarded the self-interaction terms in order to handle the singularities in the 
velocity held. They had no other choice. Except for the initial regularization step 
we took, this discarding of the self-interaction term is the primary place where the 
derivation of the equations of motion presented here diverges from the derivation 
in |YTKd9] . Discarding the self-interaction term in |YTK09] lead to contradictory 
compatibility equations for singularities of degree 2 and higher. This is one regime 
where the self-interaction terms have an impact on the dynamics in our regularized 
formulation. One of the major Endings of |YTK09j was that one could avoid these 
contradictory compatibility conditions by limiting one’s self to combinations of 
point-vortices and dipoles. Even in this limited scenario, our equations of motion 
do not match even in a regularized sense, as the vortex cores of the dipoles are 
not advected by the (singular) velocity held in |YTKn9] . Additionally, as the 
regularization parameter goes to 0 in our framework, the velocity helds become 
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singular, and the equations of motion for the F’s will explode. So we can not 
expect to observe any form of convergence to the hnite valued ODEs of |YTK09] . 


5. Conserved Quantities 


In this section we begin to touch upon some of the symplectic geometry of 
MVBs. To begin, let us consider a general vorticity distribution u. The energy is 
dehned as 


H{u) 


1 

2 


u(z)Ga{z — z')u{z')dz dz' 


1 

2 


'ijj{z)u{z)dz. 


where ^jJ = Ga * ui- The vorticity equation, Q, can be seen as an instance of 
Hamilton’s equations on a Poisson manifold. In this case the Poisson manifold 
is the space of vorticity distributions, and the Poisson bracket is the vorticity 
Poisson bracket derived in |MW83] . As the Hamiltonian is conserved by Hamilton’s 
equations, we should expect H{u) to be constant in time. Indeed, we hnd that if 
00 satishes ([^, then 


dH 

dt 


(a;) 


I {Go, * oj{z)) dtOj{z)dz 


ip{z) dtOo{z)dz 
'll) {d^oj dy'ip — dyOj dx'ip) dz 

dy - dx ) dyoodz 


By integration by parts, we can remove the partial derivates from the a;’s to find 




UJ “ 1 “ dyx 


^dz 


0 


which vanished by the equivalence of mixed partials. 

As Q is a Hamiltonian system, we can consider searching for symmetries to hnd 
other conserved quantities using Noether’s theorem. We’ve relegated the discussion 
of the relevant symplectic structure to Appendix where derivations and proofs 
of the following can be found. Here we can summarize the appendix. 

It’s simple to observe that the Hamiltonian H is translation invariant, and that 
H is rotationally invariant as long as the kernel Ga has rotational symmetry. Thus 
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we should expect there to be conserved quantities tied to these symmetries. We 
hud that the quantities 


Jlinfi:^) = 


-yu{z)dz, / xu{z)dz 


Jang(t<j) = {x + y^)u{z)dz 


are conserved. The relationship between these quantities and the symmetries of the 
system is explained in Appendix]^ Alternatively, one can observe the conservation 
of these quantities by direct calculation in the same way that conservation of energy 
was verified. 

As the MVB ansatz is consistent with ([^ we can substitute the MVB ansatz 
into the above conserved quantities, to obtain conserved quantities for the MVB 
evolution, ([IT|) and (12). We obtain the following conserved quantities: 




Jim = 2^(ri' -Tj'I,-Tj')■ 


«= Y. (- 1 )' 


Again, the first two quantities, Jang and Jun, are momenta derived from Noether’s 
theorem for the rotational and translational symmetries of the fluid. The quantity 
H is the kinetic energy of the fluid. In section we will characterize the MVB 
dynamics as Hamiltonian systems, with Hamiltonian H. 

To each individual MVB there are numerous conserved quantities which can 
be seen as a manifestation of the conservation of circulation. To show this, let 
u = ( m , v) = {dyip^ —dxip) be a time-dependent vector held which satishes ([^. The 
how of u is the diheomorphism, —)■ which sends particle labels at time 

0 to their positions at time t. If Ut is the vorticity at time t then = Uq 

is constant in time. This conservation law can be seen as a corollary of Kelvin’s 
circulation theorem |AK92j . As a consequence, the quantity 

J{t) := jut{^t{z))f{z)dz 

is constant in time for any / G By applying the change of variables 

formula and invoking the incompressibility condition, det(Zi)<I>) = 1, we hnd 


J{t) = jutiz)f{^t \z))dz. 
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This form of writing J{t) makes sense when ut is a distribution. As a result, we 
find that for a vorticity of the form (|^ the quantity 

(13) J(t)= rr(-i)“araj(/°<i>r‘)L.«.(.) 

ies 

m+n<N 


is conserved for any / G While this conservation law holds for all func¬ 

tions with compact support, /, we do not obtain infinitely many conserved quanti¬ 
ties when ujt satisfies the MVB ansatz and S is finite. This is because the expression 
on the right hand side only depends on the A^th order Taylor expansion of / at 
Zi{0) = as is illustrated by the Faa di Bruno formula. We will not 

display the Faa di Bruno formula here because it requires nearly a page of nota- 
tional definitions before to writing it down |CS96] . Nonetheless, by computing the 
cardinality of jet spaces, one would obtain card(5') independent conserved 

quantities as a result of (13). These conserved quantities can be interpreted as a 
finite dimensional manifestation of the conservation of circulation. 


6 . Moments 

In this section we present how the moments of the vorticity distribution evolve 
in time. We will find that when the vorticity distribution is that of a MVB, then 
the moments form a closed dynamical system at finite order. 

The (a, 6)**^ moment of the vorticity, u, centered around the vortex position 
e is given by 

Ff - Vifixdxdy . 


We call the integer a + h the order of the moment. For a general vorticity, the 
evolution for the nth order moments will depend on the (n-|- l)th and higher order 
moments and so we can not concoct a closed dynamical system for the moments 
of order n and less. However, this is not the case if u satisfies the MVB ansatz, 
and the points {xi,yi) are given by the locations of the jet-vortices. If u satisfies 
the MVB ansatz (|^ then 




E (-1) 


m+n 


a\b\ 


j&s 

m<a,n<b 


{a — m)\{b — n)\ 


mx,-x,r--iy,-y,Y 


ioT a + b < N and i G S'. Given the points Zi G S', one can write the moments in 
terms of the circulation strengths, the F’s. For the moment with a+b < N with 
a, 6 G N we may invert this relationship to write F™” = Fj”"’(/r), i.e. as a function 
of the moments. Invoking the motion equations for the F’s and substituting the 
relation between the F’s and the p’s yields a closed dynamical system for the p’s. 
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Remark 6.1. This relation between the T’s and the /I’s may also be important 
in the context of plasma physics, especially when one recalls that ([^ can be in¬ 
terpreted as a one-dimensional plasma model. Specihcally, phase-space moments 
of the Vlasov probability distribntion form an important dynamical link between 
Lagrangian-particle and Eulerian-continnum descriptions. The phase-space mo¬ 
ments of the Vlasov probability distribntion provide collective coordinates for 
the Hamiltonian dynamics of ensembles of particles. For more explanation of 
this property of Hamiltonian collectivization of the phase-space moments, see 
[GS9ni IHLSQOi IGHTOSbl IGHTOSa] . In plasma dynamics, the phase-space mo¬ 
ments arise from a Taylor expansion of the Vlasov particle distribntion, taken 
aronnd its centroid in phase space. For planar incompressible flow of an ideal 
flnid, the phase space comprises the (x, y) Lagrangian coordinates of a flnid par¬ 
ticle, and the corresponding moments arise from Taylor expansions around the 
centroid of the (smooth) vorticity distribution. The duality between the resulting 
spatial moments of a smooth vorticity distribution and the MVBs correspond¬ 
ing to higher-order singular vorticity distributions also obtained from a Taylor 
expansion raises the intriguing question of Ending a relation between these two 
types of dynamical description. This question is particularly intriguing because 
the dynamics of moments beyond quadratic order in general does not close to form 
a finite-dimensional Hamiltonian system, while the dynamics of MVBs closes at 
every order. 


Remark 6.2. There exist other systems for approximating the dynamics of mo¬ 
ments which differ from the one presented here. In particular, the equations of 
motion for the moments here form a closed system at order iV, whereas other meth¬ 
ods for deriving dynamical systems for moments |UWB12| INSUWO^ IGHTOSb]. 
IGHTOS^ require truncations in order to form a closed system. For example, 
[UWB12] approximates the stream function as a sum of Hermite functions with 
evolving centroids and weights. In order to obtain the evolution for the weights 
and the centroids they project the viscous vorticity equation onto this space via 
projection. The resulting formulas are explicit and efficient to compute, albeit 
more complex than the formulas found in this paper. The primary source of error 
for |UWB12] over long times is the discrepancy between the projected evolution 
equations and the true evolution equations. In contrast, we approximate an Euler 
fluid with a regularized fluid equation which we solve exactly. This is not to say 
that error is not accumulated in time. The primary source of error for our method 
over long times is the discrepancy between between the regularized fluid equations 
and the true fluid equations. 

Admittedly, the equations of motion for the moments in [UWB12] bear some re¬ 
semblance to the equations of motion for the T’s in our method. Both are quadratic 
in their respective variables, with coefficients involving combinatorial functions. A 
more precise relationship, if one exists, is difficult to discern. Philosophically, the 
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methods share much in common. However, due to the fundamental approxima¬ 
tion technique of projecting the equations of motion versus regularizing them, the 
methods are indeed distinct. This difference cascades throughout the study of 
both methods. For example, the convergence for |UWB12] is obtained via the 
convergence of spectral approximations, while the convergence of our method is a 
corollary of the convergence of a regularized fluid model (see |MM2ni,^ IFHTOlj 
for such convergence proofs). 

7. Numerical Aspects 

In this section we discuss various numerical aspects of using MVBs to model 
fluid dynamics. We will observe how MVBs can be used to reduce the number of 
necessary pairwise computations without a drastic compromise in accuracy. We 
will also present an algorithm for constructing an initial condition of MVBs from 
a given stream function. 

Remark 7.1. We refer to |CN88j for a convergence proof and error analysis of the 
1st order case. A convergence proof is beyond the scope of this article. Suffice it 
to say, such a proof would likely resemble |CN88] . 

7.1. Grouping and reduction of pairwise computations. Let us consider the 
vorticity distribution 

OJ = -|- F2(j^2 • 

If zi and Z 2 are close, we can dehne the quantities 7 = {z\\ z< 2 )l‘l and bz = z\ — z^ 
to obtain the approximation 

y u{z)f{z)dz = Tif{zi) + T2f{z2) 

= Fi (^f(z) + dJ{z)-^-^ + dyfiz)-^-fj 

+ ha (^f{z) - dj{z) ■ y - dyf{z) ■ y^ + o{h) 

where h = ||fc||. Therefore the distribution 

cD = F4 + T^dj, + rydyb, 

with 

SlJ 

F = Fi + F 2 , F" = y(F2-F0 , F^ = ^(F 2 -Fi) 

serves as a o{h) approximation of uj in the sense of distributions. Moreover, the 
stream fnnction -0 := * a; is an o{h) approximation of "0 := Gs * oo in the 

traditional sense of analysis on functions. 
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We have just described the hrst case of grouping two iVth order MVBs concen¬ 
trated at Zi and Z 2 into a single {N + l)th order MVB concentrated at the average 
position More generally, we can consider the ansatz 

a;= ^ 

m+n<N 

and observe 

jw(z)f(z)dz^ Y, (-ir+"(rrarapLi) + rr9"apL2)) 

m+n<N 

Y (-L^+’T’r (srs;/(s)+sr‘ap(L)) ■ y+srsr‘/(L)) • y 

m+n<N 

+ (-ir+"rr (d':d'‘f{z) - ar'spd.-)) • y - src'AL)) ■ f) | 

+ o{h). 

The above computation implies that the quantity 



u : = 


E T" 


m-\-n<N -\-l 


I pmn 

“T -L 2 


Snr' 


m—l,n 


1,71' 


/pm,n—1 


-r 


m,n—1' 


nr- 




serves as an o{h) approximation of uj. Of course, this again implies that the 
corresponding stream functions are approximated to order h as well. Note that uj 
is concentrated above a single point, z, while oj is concentrated above two points. 


Remark 7.2. Such reductions are even more dramatic when considering higher 
order jets. In particular, 2^ zeroth order MVBs can be approximated with a 
single iVth order MVB by applying the above approximations iteratively. 

The computation of pairwise interactions in the vortex method was once a ma¬ 
jor bottleneck in implementing the standard vortex method for real-world applica¬ 
tions. It was not until the invention of the fast multipole method, that it became 
tractable to compute millions of pairwise interactions by reducing the complexity 
from an 0{'n?) calculation to an Cl(nlog(n)) calculation, where n is the number of 
vortices |GR87] . However, in the case of viscous fluids with boundaries, vorticity is 
shed from the boundaries. As a result, the vortex blob method of |Cho73] created 
new vortices at the boundary by using the Kutta condition as a creation criteria. 
For these applications, n will grow in time without bound, and some means of 
discarding vortices must be invoked. It is here that the grouping of MVBs could 
be useful. If one merges two Vth order MVBs to obtained a {N -|- l)th order MVB, 
the amount of scalars and data typically increases. So one must still make a tough 
decision as to what data to discard (e.g. through some tolerance or by simply 
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truncating at level M). Nonetheless, the analysis presented here could shed light 
on how best to implement this approach. 


Remark 7.3. The merging of blobs of vorticity has been studied analytically |MZM88] 
and numerically |WM93irMZM88ilLDV02j . as well as in the laboratory |FDMM9T| . 
All of this study has been in the slightly viscous (or nearly inviscid) regime. The 
grouping approach discussed here can be used to numerically resolve such collision 
events. In theory, there is no issue with collisions because we are considering regu¬ 
larized vortices where the induced velocity field from a single MVB is always finite. 
However, as 5 becomes smaller, the velocity near the vortex core diverges. This 
should be of concern as the convergence analysis of the vortex method pre-supposes 
that (5 1. Typically such a near collision is handled by using a smaller time-step 

(as the ODE is quite stiff). Grouping of MVBs suggests an alternative by avoiding 
this pair-wise interaction altogether. Perhaps such an approach could be viewed 
as a variation of the punctuated dissipation events described in |WM93j where 
an initial vorticity distribution is found to asymptotically approach a smoother 
axisymmetric vortex blob, and discrete vortex mergers are implemented to model 
this behavior. 


Remark 7.4. There are qualitative questions which arise from mergers. For exam¬ 
ple, when two 0th order vortex blobs are near each other, they will typically scatter 
after some hnite time. Merging these blobs into a single 1st order blob will prohibit 
this scattering event from ever occurring. That both the zeroth-order MVB solu¬ 
tion and the merged 1st order MVB represent exact solutions of the fluid (after the 
merger event) is attributable to the long term sensitivity to initial conditions near 
collision events. The scattering angle can be virtually anything since zeroth-order 
MVBs can waltz around each other many times before scattering. The amount of 
time two zeroth order MVBs can spend waltzing around each other, and perhaps 
the merged solution represent some sort of limiting solution. That is to say, the 
merged solutions can be interpreted as the “waltzing for eternity” solution. 

The irreversibility of merging is disturbing when one takes it to its extreme, one 
massive high order MVB. In order to address this, a means of splitting high order 
MVBs into lower order ones should be considered. The primary difficulty here is in 
determining when to split. In the case of mergers, we can decide to merge MVBs 
when they are close. Such a criterion is not immediately apparent in the case of 
splitting MVBs. 


7.1.1. A numerical experiment with grouping. For illustrative purposes we can nu¬ 
merically group four 0th order MVBs into two 1st order MVBs, and then one 2nd 
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order MVB. In particular, we can consider the initial condition 


Zo = 

(-0.25,-0.25) , 

ro = 

0.3 

Zi = 

(-0.25, 0.25) , 

ri = 

- 0.35 

Z2 = 

( 0.25, 0.25) , 

F2 = 

- 0.2 

Z 3 = 

( 0.25,-0.25) , 

r3 = 

0.4 


The corresponding dynamics are depicted in the top row of hgure 

Next we group Zi with zq and Z 2 with z^ in order to obtain two 1st order MVBs 
with initial condition 


(15) 


zo = (-0.25,0.0) , To = -0.05 
zi = ( 0.25,0.0) , Ti = 0.20 


= 0.0 , rs' = 0.1625 

= 0.0 , rs' = 0.15 


The corresponding dynamics are depicted in the middle row of Figure The 
dynamics appear qualitatively similar at the beginning of the evolution. Then the 
dynamics diverge around time t = 150 when the two 1st order MVBs separate 
from one another, in contrast to the dynamics of the 0th order MVBs. 

Finally, we group the two 1st order MVBs to obtain a single 2nd order MVB. 
Again, the dynamics appear qualitatively similar at the beginning of the the evo¬ 
lution. Oddly, the dynamics of the 2nd order MVB appear qualitatively similar to 
the 0th order case even at t = 253. As there is only a single vortex, the separa¬ 
tion of vortices mentioned in the 1st order MVB experiment is not possible here. 
As a result the dynamics of the original 0th order MVB dynamics appears to be 
approximately recovered. 


7.2. Approximation of initial conditions. In this section we will illustrate 
how initialize MVBs when given a stream function ^ at time 0. We can begin by 
dehning an inner-product on the space of distributions on given by 

{uJi,UJ2)gs ■= j 0 Ji{z)G5{z - z)uj2{z)dzdz. 

Consequently, the energy of the fluid is given by H{u) = = \{u,u)gs- 

Let A be a compact set and let 0 < h 1 be small so that we may dehne 
the hnite grid A/^ = {{ah^hh) G A | (a, 6) G Z^}. Q Given an a; G 'D'(M^), we 
can attempt to approximate oj via Dirac-deltas supported on hi?. There is a 
natural way to do this with respect to the inner product {■,-)gs- We could dehne 
by requiring the error, cn® —ca, to be (•, •)G'a-orthogonal to 6^ for 
each 2 ; G A/^. This means that Gs * oj{z) = ^i^siz — Zt) for each z G A^. Thus 
— Zi) can be seen as a 0th order approximation to ?) = G^^u 

^ The choice of K should depend on the initial circulation w, e.g. if the u has compact support 
than any K which contains the support of oj would be a good candidate. Nonetheless, having to 
choose AT is a weakness of the given approximation procedure. 
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Figure 4. Snapshots of the evolution for various initial conditions 
at time t = 0, 51,101,152, 202, 253. The top row depicts the evolu¬ 
tion of four 0th order MVBs given by the initial condition (14). The 
middle row depicts the evolution of two 1st order MVBs obtained 
by grouping. The bottom row depicts the evolution of one 2nd order 
MVB obtained by grouping. 


because = ipiz) for all 2 ; G Ah. Therefore, for smooth oi’s, we obtain an 

error of order 0{Ax) for a grid-spacing of Ax using 0th order MVBs. 

The same reasoning applies if we consider = '^im+n<N^T^^T9ySzi. We 
dehne the scalars T™"’ via the equations 

- zj) 

j 

for = Gs*uj, Zi e Ah, and |/3| < k. Then 'il)h\z) = 

Zi) serves as an order k approximation of "0 when -0 G C^. In particular, for smooth 
oi’s, we obtain an error of order for a grid-spacing of Ax using /cth order 

MVBs. 

As an example, we numerically compute the corresponding approximations of 
the stream function 

(16) il){x,y) = exp(—r^) — exp(—r^/2) 

The results are depicted in Figure where we observe sup-norm convergence on 
the interior of K. In particular, we measure the sup-norm error on the subregion 
(—3 < x,y < 3) with K = {{x,y) \ —6<x,y< 6}. We observe convergence 
using MVBs at orders zero, one, and two. In each case, a grid spacing is reached 
where the error plateaus (possibly due to machine precision). Nonetheless, higher 
order MVBs appear to out perform lower order ones for smaller grid spacings. In 
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particular, we observe slopes in a log-log plot of magnitudes 1,2, and 3, suggesting 
that 1st, 2nd, and 3rd order convergence rates for 0th, 1st, and 2nd order MVBs 
respectively. 



Figure 5. A convergence plot of the error in the sup-norm of the 
reconstructed stream function approximated using MVBs of order 
0 , 1, and 2. 

In terms of complexity, in order to achieve a desired error bound, etoi > 0, you 
would need to use a grid with MVBs. While the number of MVBs 

drops as k increases, one could object that a high-order MVB is much more complex 
than a low order one. However, the number of degrees of freedom for a kth order 
MVB is 2 -|- which monotonically converges to a constant (roughly 

9.39) as /c —)■ cx). Therefore the number of degrees of freedom is dominated by 
0 (e,as well. In other words, when -ip is highly differentiable we observe 
benehts in terms of complexity and storage to using a larger k regardless of weather 
one measures complexity by the number of parameters to keep track of, or the 
number of MVBs. 


8. Numerical experiments 

In these section we present the results of numerical experiments involving small 
numbers of vortices, for V = 1, 2, and 3. 

8.1. Behavior of isolated MVBs. Next, we will briefly explore the behavior of 
a single isolated kih order MVB with T™'"' = 0 with m + n < k for k = 0,1, 2. This 
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case allows us to investigate the dynamics induced by the higher order circulation 
variables in the absence of the lower order ones. 

8.2. Order 0. The behavior of a single 0th order MVB is explicitly solvable be¬ 
cause the dynamics are stationary. 

8.3. Order 1. The behavior of a single 1st order MVB with T = 0 is explicitly 
solvable. Given the initial condition (x(0), i/(0), r(0), r^(0), r^(0)) with r(0) = 0 
we hnd 

x{t) = x{0) + v^t, y{t) = y{0) + v^t, r(f) = r(0) 
r^(f) = r^(o), Ty{t) = T\o) 

where = r"'( 0 )da;,/G 5 ( 0 )-hr 2 ^( 0 )^yG 5 ( 0 ) anduy = -r 2 '( 0 )^G 5 ( 0 )-r*( 0 )dj, 2 ;G 5 ( 0 ). 
In Figure we depict such a trajectory with initial condition 

(17) a;(0) = -3, i/(0) = -3, r(0) = 0, T" = 1, T^ = 1 



Figure 6. A 1st order MVB with initial conditions given by 0 
with snapshots taking at f = 0,10, 25. 


8.4. Order 2. The behavior of a second order vortex does not seem to be explicitly 
solvable. Here we consider initial conditions for which 


(18) x(0) = 0, 1/(0) = 0, F(0)"" = 1 

and all the other circulation variables are initially set to 0. The results are depicted 
in Figure We observe a structure which rigidly rotates counter-clockwise. 


8.5. A scattering expiriment. Next we consider two MVBs. The first is a first 
order MVB with an initial velocity pointed just slightly above origin. The second 
MVB is a standard zeroth order vortex located at the origin. Specificaly, we 
consider the initial conditions 


zo = (20.0, 0.25) 
= ( 0.0,-0.25) 


F°’° = 0.0 , Fo’° = 0.0 , F°’^ = -1.0 

Fo’° = 1.0 , Fo’° = 0.0 , Fq’^ = 0.0 
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Figure 7. A 1st order MVB with initial conditions given by (18) 
with snapshots taking at t = 0, 5,10,15,20, 25. 


with r™" = 0 for m + n > 1 and i = 0,1. The vortex at the origin appears 
to remain at the origin throughout the numerical run (t = 0 to t = 150). The 
hrst order MVB starts by moving to the left in a straightline until it comes into 
proximity of the zeroth order vortex. Then the hrst order MVB swings around the 
the zeroth order vortex, traversing an angle of roughly 30 degrees before zooming 
off into the lower left quadrant of the plane in a straight line. These results are 
depicted in Figure 



Figure 8 . A numerical run is shown with mirror image initial 
conditions for two 1st order MVBs, as given in (19). From left 
to right and top to bottom these are snapshots at times t = 
0, 25, 50, 75,100,125 respectively. 


8.6. The method of images. Here we incorporate hrst order MVBs into the 
method of images jJac75( ILSll] . We consider the initial conditions consisting of 
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two first order MVBs which are mirror images of each other with respect to the 
x-axis. By symmetry, the resulting vector-held should be tangential to the x-axis, 
and provides a means of considering a boundary that satished the no-penetration 
condition. Specihcally, we consider the initial condition: 

1^0 = (1.5, 1.5) , rr= 0.5 , ry= 0.5 , rS’‘ = i.5 

Pi = (i.5,-1.5) , rJ“ = -0.5 , ry = -0.5 , rw^i.s 

with r™" = 0 for m -|- n > 1 and i = 0,1. 

The resulting dynamics depicted in Figure shows that as a hrst order MVB 
approaches a boundary it will turn its motion along the boundary and then move 
away so that its angle of rehection equals its angle of incidence. 



Figure 9. Numerical results are shown for two hrst order MVBs 
with mirror-image initial conditions given by (20). From left 
to right and top to bottom these are snapshots at times t = 
0,1.7, 3.4, 5,6.7, 8.4 respectively. Apparently, a hrst order MVB re- 
hects elastically from a hxed boundary, so that its angle of rehection 
equals its angle of incidence. 


9. Hamiltonians and symplectic structures 

In modern Hamiltonian mechanics, as described in |AM78( lArnOO] , the Hamil¬ 
tonian is a function on a symplectic manifold, which produces equations of motion. 
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An important instance of a symplectic manifold is a coadjoint orbit (defined be¬ 
low). In this section we compute the coadjoint orbit of a MVB as well as the 
associated symplectic structure. The coadjoint orbit of an initial vorticity distri¬ 
bution ojq comprises the set 

Orb(a;o) := {wq o \ (p G SDiff(M^)}. 

In fact Orb(ci;o) inherits the structure of a smooth manifold, and a tangent vec¬ 
tor on Orb(ci;o) at the point Cj G Orb(ci;) is given by a distribution of the form 
T^[a)] := udxOJ -f vdyUJ for some (non-unique) divergence free vector held u = 
iu^v) G Xdiv(lR^)- The symplectic structure is nothing more than a special case of 
the one derived via the Kirillov-Kostant-Souriau theorem |AM781 see the boxed 
formula on p.303]. In particular, the symplectic structure on Orb(a;) is given by 

( 21 ) T^Jca]) = j u{z){ui{z)v2{z) - Vi{z)u2{z))dz. 

When a; is a smooth distribution, the symplectic structure may be identihed with 
a differential 2-form and this formula matches the symplectic form derived on page 
313 of |MW83j . In the case that oj satishes the ansatz ([^, we hnd that given any 
ip G SDiff(M^) that 

j UQ{ip~^{z))f{z)dz = j uo{z),f{ip{z))dz 
= 7Tda\z=zAf 

Here we have used the change of variables formula and the fact that det(Zl(p) = 1. 
By the multivariate Faa di Bruno formula, the expression da\z=Zi{f ° ^){^) is a 
sum of the partial derivatives of / at the points ip{Zi) of order less than that of 
the multi-index a |CS96j . Thus uq o (p~^ is contained in the hnitely parametrized 
subset := for any ip G SDiff(M^). Therefore Orb(a;o) is a 

hnite-dimensional manifold when cuo satishes the jet-vortex ansatz. 

Having identihed a symplectic manifold, Orb(a;o), we can then ask the question 
“are the dynamics Hamiltonian on Orb(a;o)?” Of course, the answer is “yes”. 
This is the primary content of |MW83j . We provide our own explanation here for 
convenience. 

For a general vorticity distribution ca, we may consider the kinetic energy Hamil¬ 
tonian 

( 22 ) H{oj) ~ 2 J~ 

Where oo may be of the form (|^. In order to hnd Hamilton’s equations on Orb(a;o) 
choose some oj G Orb(a;o) and calculate the vector Xh{oj) tangent to Orb(a;o) 
given by Hamilton’s equations. It must be the case that Xh{oj) = £ii[oj] for 
some (non-unique) vector-held u = (m,u) G Xdiv(lR^)- Our goal is to solve for 
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u. By the definition of the Hamiltonian vector field Xh we see that for any 
u' = {u',v') e Xdiv(IR^) 

/ /* 

u(z) (u(z)v'(z) - v(z)u'(z)) dz = £ii'[u]) = - 

= — J Gs{z — z)u{z) {£^i[ijj\) {z)dzdz. 

If we let -0 := * a; = / G 5 (- — z)uj{z)dz then integration by parts implies 

j u{z)£it'['i/j]{z) = J u{z) {u\z)^x^/J{z)+v'{z)dyi>{z))dz. 

We see that u = {—dy'ip, is one possible solution. As H is non-degenerate 
on the tangent spaces of Orb(c(;), this is the unique solution. As a result, the 
evolution prescribed by is precisely Q. This proves that Q can be seen as 
a Hamiltonian equation on Orb(a;o) with respect to the symplectic structure (21) 
and the Hamiltonian (22). 


9.1. The first order case. Let us illustrate these Hamiltonian results for the 
case of the first order MVB. Let 2 ;i,..., G be distinct and define the initial 
vorticity distribution 


N 

i=l 

We desire the to determine the coadjoint orbit, Orb(a;o), and the symplectic struc¬ 
ture. 

Indeed, we find that for any function / 


Juj{ip \z)) f{z)dz := juj{z)f {ip{z))dz 

- ltdx^''\z=zAf\z=^{Zi) - 7!da:(py\z=zAf\z=V^iZ,) 

- lldy£%=zAf\z=^{z,) - lldy^y\,=zAf\z=^{z^) 

Collecting like terms we find 

ujoip ^ + £idxdy,(z,) + £^idy6^(^Zi) 


where 


r = 



D(p{Zi) ■ 
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By varying ip we can obtain any collection of distinct points zi,... ,Zn G and 
any collection of non-zero vectors Fi,..., G This proves 


Orb(i^„) = \ + rf4*, + rfdyi,, I z, e K=, (rf.r?) e K^{o} 


2=1 


= {(zi,... ,z„,Fi,... ,F„) I Zi e M^Fi e ]R^\{0}, (i 7^ j Zi^Zj)}. 


To derive the symplectic structnre recall the symplectic structnre for a general 
vorticity (21). Now let u = + TfdxS^^ 


+ rpySy 


In this case the left hand 


side of (21) can be compnted with respect to divergence free vector field u = (n, v) 
and u' = iu' v') as 


j Uj{z) {uv' - vu') {z)dz = '^i{u{Zi)v'{Zi) - v{Zi)u'{Zi)) 

T T uv U ^xU U V^x'}\z=zi 

-h T^iiu^yv' + UV\y - u\yV - u'v^y)\z = Zi 


Note that this is written entirely in terms of the 1st order Taylor expansion of 
u and u' evalnated at Zt. Moreover, £ii[uj\ = '^ou{zi)dx5zi + ... also has the 
property that it only depends on the first order Taylor expansion of u and v at the 
points , Zn. Therefore, both sides of (21) can be written as a fnnction of the 

hnite collection of numbers u^Zi)., Du{zi),v{zi), Dv{zi). The result then follows by 
identifying the scalars 


u{Zi) ^ 

u%{z.)Tt + ul{z,)Ty^tt 
uy^{z,)T^ + uyy{z,)T^, ^ ti 

This proves that the symplectic structure on Orb(ci;o) is more concretely written 
as 

n{{z,t),{Sz,5T))='yi{xi-6yi-5xi-yi) 

■ Sy, - ty ■ Sx, + dry • x^ - ^Ff ■ y, 

In essence, we have determined a finite-dimensional Hamiltonian system whos 
solutions solve ([^ when ip is related to u via an appropriate regularization. 

Remark 9.1. The use of this symplectic structure shows that the map (^j, Fj, Ff, F^) i—)■ 
u G Orb(a;) is a symplectic momentum map. 


Remark 9.2. The corresponding Poisson bracket can be represented in tabular 
form by: 
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X 

y 

r 


Yv 

X 

0 

1 

0 

0 

1 

y 

-1 

0 

0 

-1 

0 

r 

0 

0 

0 

0 

0 


0 

-1 

0 

0 

1 


1 

0 

0 

-1 

0 


The way to use this table is as follows. Let H = H{^) be our Hamiltonian 
where ^ = (a:i,a:„, 2 / 1 , Ti,r„, T^,T^, T^,T^). Hamilton’s 

equations are then given by given by 


dt ^ ’ 

hj 


for any function /, where denotes the corresponding entry of the table. In 
particular, when / = one recovers the equations of motion for the dynamics of 
the positions and strengths for a set of n 1st order MVBs. Poisson geometers call 
B'’^ a Poisson tensor |AM78] . 


10. Conclusion 

In this paper we have considered a generalization of the standard vortex blob 
method, obtained by augmenting the vortices with higher order circulation vari¬ 
ables and dubbing them multipole vortex blobs (MVBs). By viewing the vorticity 
equation as an advection equation, we have obtained equations of motion for these 
MVBs. 

The extra degrees of freedom of MVBs resulted in richer dynamics near the 
vortex core. Moreover, these new vorticity carrying elements exhibited a variety 
of novel types of solution behavior. We also observed faster convergence rates in 
space using higher order MVBs. Moreover, we proposed a scheme to decrease the 
number of pairwise interactions, by grouping MVBs of lower order into a smaller 
number of MVBs of higher order. Lastly, the implications of Kelvin’s circulation 
theorem were substantially richer in the case of MVBs than they were for the 
standard vortex blob method. 

We have demonstrated the behavior of the MVBs with a sequence of simple 
numerical experiments consisting of small numbers of MVBs of various degrees. 
We found that 1st order MVBs correspond to sums of vortex blobs and regularized 
dipoles which simply propagate themselves forward, while the 2nd order circulation 
variables activate richer (non-propagating) dynamics near the vortex core. 

Finally, we derived the symplectic structure of MVBs using methods from |MW83j . 
The resulting structure turned out to be a direct generalization of the standard 
symplectic structure for vortex blobs. 

The multiscale nature of ideal fluids is the principal obstacle to obtaining ac¬ 
curate models |Cho94( Ch. 3]. The use of MVBs augments the standard vortex 
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blob method by allowing for singular vorticity distributions which model dynamics 
below the regularization length scale (i.e. at order 6^^ with 5^1 for a kth order 
jet-vortex). As the dynamics of MVBs are relatively easy to derive, and their 
analysis is tractable, we believe that MVBs will be of considerable value in un¬ 
derstanding the place of regularized fluid models within the computational fluids 
community at large and they should provide renewed interest in the vortex blob 
method. 

Future avenues of inquiry could include: 

• MVBs on manifolds, such as the sphere 

• The convergence properties of the MVB method 

• How does one choose the regularization length-scale in relation to the grid 
resolution. This relationship is addressed quite well for zeroth order MVBs 
in jBM82j . It is not clear if higher order MVBs change those results. 

• An investigation of the kinetic theory of MVBs. 
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Appendix A. Distributions 

The vorticity, cu, should be viewed as a distribution and the term '^dxOj” should 
be viewed as a distributional derivative. When u G is a smooth distribution 

there is little harm in naively interpreting a; as a smooth function on However, 
when uj is not smooth (e.g. a Dirac delta distribution), then one needs to invoke 
the mathematics of distributions as distinct from that of real valued functions. 
Therefore, we have included this appendix to remind the reader of the basic theory 
of distributions. The main reference for this section is [HorOd] . 

The space of distributions P'(M^) is the dual-vector space to the space of smooth 
functions with compact supper (^“(M^). Therefore a distribution is dehned by how 
it maps functions to real numbers. 

The distributional derivative of a given distribution u G P'(M^) in the Ah coor¬ 
dinate direction may be dehned as the distribution diU obtained by 

J diUj{z)f{z)dz = - J uj{z)dif{z)dz. 

For example, the Dirac-delta distribution, (5o, is dehned as the unique distribution 
such that 

j Mz), f(z)dz = /(O), V/ e 
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The distributional derivative, di5o, is given by 

j dAiz)fiz)dz = -d,f{0) , V/ e Co“(m2). 

Given a distribution u G D(M^) and a function g G one can define the 

distribution gu as 

i 9 (^)iz)fiz)dz = Ju{z)g{z)f{z)dz , V/ G 

For example, we hnd that g6o = g{0)6o. A slightly more involved, but standard, 
example is given by the computation of g dido. We hnd 

j {gdi6o){z)f{z)dz = j di6o{z)g{z)f{z)dz = -g{0)dif{0) - dig{0)f{0). 

Therefore 

gdido = g{0)di6o - dig{0)6o. 

On the left hand side, note that ^'(O) and dig{0) are merely real numbers, which 
are multiplying the distributions dido and 5o. More generally, we hnd 

j ( 3 d^d;s„){z)f{z)dz = (-i)’"+"ca^(/g)(0) 

= (-1)“ ^ (“) (") (c^9p(o)) (arvr‘9(o)) 

£^k=0 ^ / V / 

which means 

m,n 

( 24 ) gdTd;;do = 

e,k=o 

Appendix B. Symmetries and conservation laws 

The main reference for the material presented in this section is |AM78] . Let 
G be a Lie group with Lie algebra 0 . We will denote the dual of 0 by 0 *. A 
(left) group action of G on a manifold S' is a map p : G x S' —)■ S' such that 
g{gh, x) = p(p, p(h, x)) for all g,h,^ G and a; G S. 

Remark B.l. The group action q is not to be confused with the huid density, often 
denoted as p in huid mechanics. This appendix relates to more general mathe¬ 
matical constructions which are useful, but not necessarily in the usual purview of 
huid mechanics. In particular, the symbol q is the Greek letter ‘r’, which refers to 
the word “representation” as in “representation theory”. 

One can also construct a group action, Dg : G xTS ^ TS, given by Dg{g, u) := 
^ 1^=0 ~ ^ ^ 9 ^ G. There is also a natural Lie-algebra 

action, which one could also denote hy g : q x S ^ TS dehned by g{^,x) := 
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^|^_o Q{ge, x). In particular, the map •): S' —?• TS" is a vector field on S which 
we call the infinitesimal generator of When no confusion arises, it is typical 
to use the notation g ■ x, g ■ u, and ^ • x to denote g{g,x),Dg{g,u), and g{f,x), 
respectively. Finally, if (S, fl) is a symplectic manifold, then we say that G acts 
symplectically (or canonically) if Qg.x{g-u, g-v) = Qx{u,v) for all g E G, u,v E T^S 
and X E S. 

Let us now recall the notion of a momentum map |AM78( Definition 4.2.1]. 
Given a symplectic manifold (S', D) and a Lie group G which acts on (S', D) sym¬ 
plectically, a momentum map is a map J : S —?■ g* such that 


(^{3, f) — i^s^ , 

where (J,^ denotes the real-valued function on S' obtained by pairing J with an 
arbitrary ^ G g, and where G X{S) is the infinitesimal generator of ^ on S'. 
Equivalently, we could express the previous condition as 


(25) 


d(J,.:))(x) ■ Vx = ■ x,v^) 


for all ^ G g, X G S' and E T^S. Momentum maps are significant for a number of 
reasons. In particular, given a Hamiltonian on S' with G-symmetry, the momentum 
map J will be conserved under the evolution of Hamilton’s equations |AM781 
Theorem 4.2.2]. This is the Hamiltonian version of Noether’s theorem. 

In our case S = Orb(a;o) = {tuo o (p~^ \ ip e SDiff(M^)} is a coadjoint orbit of 
some vorticity distribution on Tangent vectors on S' are of the form 

Xxfoj] := udxOJ + vdyOJ 

for a (perhaps non-unique) divergence free vector-field u = {u,v). Under this 
identification, the symplectic form at some cu G S is given by an application of 
Kostant’s formula |AM78] . This is derived in section]^ and found to be 

D^(T^Jw], T^Jw]) := / uj{z){ui{z) X U 2 {z))dz. 


where Ui and U 2 are divergence free vector-fields and x denotes the planar cross 
product. Here we interpret the planar cross product as taking values in the space 
of real-numbers so that Ui x U 2 is merely a smooth function. 


We now will translate formula (25) to this more specific scenario. Assume G 
acts upon then G also acts upon distributions and upon S by symplectic group 
actions. In this context, a momentum map J associated to a G-action is defined 
by the equation 


(26) 


yj.o 

5ijj 


z) {£^[uf\) {z)dz = / uj{z) ((^ • z) X u(z)) dz 


Where ^ ■ z denotes the action of G g on x G and “x” denotes the planar 
cross-product. 
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B.l. Translational symmetry and Jun. The group acts upon by transla¬ 
tion. That is to say, by sending any 5 G to -f 5 G for any 5 G This fact 
induces an action on smooth fnnctions. In particular, there is a natural (right) ac¬ 
tion on sending the function cj){z) to the function {z)*(f){z) := (j){z — z). We 

denote the inverse operation by {z)^(j){z) := cl){z + z). This induces a (left) action 
on distributions which sends u G 'D'(M^) to the distribution {z)^uj{z) := uj{z + z). 
As a translation of by 5 is a volnme-preserving diffeomorphism, we see that the 
coadjoint orbit S C 'D(]R^) is invariant nnder this action. Moreover, we observe 
the action, restricted to S, is symplectic because 

(-^^[^] 7 [^]) 7 [^*^]) 

= J z^u{z) {z^{u X v)) {z)dz 
= J u}{z + z) {u X v) {z + z)dz 
= j uj{z) {u X v) {z)dz 


Given this symplectic action, we can seek a momentum map, Jun : S —)■ (M^)*. 
Consider an arbitrary element of the Lie-algebra 5z = {6x,6y) G and use (26) 
to obtain 

f 6{Jlin,Sz) 


6u 


{z)£it[u]dz = / u{z) {6xv{z) — 6yu{z)) dz 


We can re-write the right hand side as 


and npon integrating by parts this is equivalent to 


J £j 2 [uj]{z) {5xy — Sy x) dz 


Therefore, “cancelling” the arbitrary vector £v[u:] from both sides we hnd 


Integrating by oj we find 


Soj 


6yx — 6xy 


Jiin(t^) = - / uj{z)ydz, 


u 


{z)xdzj 


d2^* — 


If uj satisfies the MVB ansatz 


a; = + T^dyd, 


+ 
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then 

Ju,(a;) = 5^(r^-r,2/„r,a:,-rn. 

i 

The terms of the MVBs beyond the first order do not inflnence Jun- 


B.2. Rotational symmetry and Jang- The gronp SO(2) acts npon by rota¬ 
tions abont the origin sending z & ^ to Re ■ z ■.= {cos{9)x — sm{9)y,sm{9)x + 
cos{9)y). For any 9 G SO(2), there is a natnral action on sending the fnnc- 

tion cj) e to the fnnction 9*(j){x, y) := (/)(cos( 6 ')x—sin( 6 ')j/, sin( 6 *)a;-l-cos( 6 ') 2 /). 

The corresponding action on vector-helds and all other objects on follows natn- 
rally. In particnlar, the left-action on distribntions sends u G to the distri- 

bntion 9^uj{z) := u{Re-z). Again, we can verify that the coadjoint orbit S is invari¬ 
ant nnder this action, and that SO(2) acts symplectically upon S through compu¬ 
tations which are analogous to those performed in the previous subsection. Given 
this symplectic action, we can seek a momentum map, Jang : S —?• so(2)* = M. 
Consider an arbitrary element of the Lie-algebra ^ G so(2) = M and use ( [26| ) to 
obtain 

( ang( — f {—y^V — X^U) dz 

duj J 



We can re-write the right hand side and integrate by parts to hnd 


= -C / 


-I- y'^ 


dz = ^ J £4^]{z) ^ 


-I- 


dz. 


As the vector £it[uj\ G T^S is arbitrary we find 


^(Jang,0 .x‘^ + y‘^ 


*ang 

5ljJ 


= i- 


Integrating by u we hnd 


'angv-/ 2 

If uj satishes the MVB ansatz 


Jang(l^) = - uj{z){x^ + yjdz G M = so(2) 


a; = + VWJz. + R... 


then 


>ang 


+ vl) - rfi. - rp. + rr + r 


yy 

i ■ 


The terms of the MVBs beyond the second order do not inhuence the angular 
momentum Jang- 
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